Non-equilibrium Phase-Ordering with a Global Conservation Law 
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In all dimensions, infinite-range Kawasaki spin exchange in a quenched Ising model leads to an 
asymptotic length-scale L ~ (pt) 1 ^ 2 ~ t 1 ^ 3 at T = because the kinetic coefficient is renormalized 
by the broken-bond density, p ~ L~ x . For T > 0, activated kinetics recovers the standard asymptotic 
growth-law, L ~ i 1//2 . However, at all temperatures, infinite-range energy-transport is allowed by 
the spin-exchange dynamics. A better implementation of global conservation, the microcanonical 
Creutz algorithm, is well behaved and exhibits the standard non-conserved growth law, L ~ i 1 ' 2 , 
at all temperatures. 



Globally conserved dynamics offers a useful test of our 
theoretical understanding of phase-ordering jlj because 
it allows access to off-critical quenches without requir- 
ing diffusive "modcl-B" transport of the order-parameter. 
RG ||] arguments show that adding global conservation 
changes neither the critical dynamics nor the asymptotic 
growth law after a quench into the ordered phase. This 
has been confirmed by numerical studies Q. The only 
change to non-conserved dynamics is to impose dt{4>) = 
where, using spin language, (f> is the local magnetiza- 
tion (((j)} is its spatial average). The standard "model- 
A" dissipative dynamics are dt<fi = —TSF/Scf), where 
F[{(j)}] = J d d r[{\7(j)) 2 + O 2 - I) 2 ]. The global constraint 
is imposed with a a uniform magnetic field term, H(<f>), 
added to F. The field H(t), effectively a lagrange multi- 
plier, is decreased smoothly while the ordering progresses 
so as to maintain (</>) constant. 

With a global conservation law, ((f)) is a tunable pa- 
rameter that does not change the growth-law but does 
affect scaled correlation functions and other aspects of 
the system. One can use this to explore what "universal- 
ity" entails in non-equilibrium phase-ordering systems. 
For example, with global conservation in Ising models, 
Sire and Majumdar jlj have shown that the autocorrela- 
tion exponent A depends on the net magnetization. This 
author [EJ has found spatial anisotropy in scalar phase- 
ordering systems that depends strongly on (</>). 

However, there has been uncertainty about how to im- 
plement global conservation. Tuning a magnetic field to 
maintain (<f>) is impractical in computer simulations, due 
to the stochastic nature of most algorithms. Two al- 
ternative algorithms have been used for Ising models: 
a Creutz demon and infinite-range Kawasaki spin- 
exchange ■ Differences between these two imple- 
mentations at low temperatures has led to some confusion 
in the past In this report we clarify the situation. 

For a 2d Ising model on a square lattice, consider 
infinite-range Kawasaki exchange dynamics HI]]?]]: two 
randomly selected spins are exchanged under heat-bath 
dynamics, i.e. with probability 1/[1 + cxp(AE / ksT)], 
where AE is the energy change under the spin exchange. 
The exchange satisfies global conservation. However 
at low temperature the asymptotic dynamics are acti- 



vated, and cross-over to a different growth law at zero- 
temperature. 

At T = 0, no exchanges that increase the energy are 
accepted. Starting from random initial conditions, aggre- 
gation of spins takes place and the broken-bond density 
p(t) will be proportional to the domain wall density at 
late times, p ~ L~ x . Isolated spins with all bonds broken, 
of number-density N, will aggregate rapidly onto exist- 
ing domains, dtN cx —pN, and may be ignored. Spins 
which are part of domain walls can exchange, but only 
with the fraction p of the system that has broken bonds 
(i.e. that are also on domain walls). This rescales the 
effective kinetic coefficient by a factor of p ~ L _1 . For 
T = infinite-range Kawasaki exchange we thus expect 
L ~ {pt) 1/2 ~ (t/L) 1 / 2 ~ t 1 / 3 , where t 1/2 is the standard 
non-conserved growth jlj . This applies in any dimension. 

For T > 0, thermal fluctuations provide partners 
for spin exchange and also allow exchange with bulk 
spins. These activated processes will dominate after the 
broken-bond density becomes comparable to the equi- 
librium average. After a quench, this will result in a 
crossover from intermediate L ~ i 1 / 3 growth to asymp- 
totic L ~ i 1 / 2 growth. [The length-scale is extracted 
from the scaling of the spherically-averaged correlation 
function C(r,t) = (<p(x)4>(x + r)) = M 2 f(r/L), where 
M is the equilibrium bulk magnetization. L is chosen 
so that /(l) = 1/2.] Results for quenches to T c will be 
similarly affected — the dynamical exponent will be un- 
derestimated while p relaxes towards equilibrium || . 

A renormalized kinetic coefficient alone will not affect 
scaled correlations. However Kawasaki exchange also al- 
lows infinite-range energy-transfer through the system. 
For example, at T = 0, an isolated spin (four broken 
bonds) may hop freely through the system by spin ex- 
change with bulk spins (no broken bonds). Similarly, a 
"bump" on a flat interface (with one broken bond) may 
hop to any other flat interface of the system. Effects of 
this energy transport might show up at low-temperatures 
in non-universal effects such as spatial anisotropy in the 
late-stages of the quench || . It is unclear if this would 
affect the autocorrelation results of Sire and Majumdar 
0, who used this algorithm at T = 0. 

A faster and safer algorithm was proposed by Creutz 
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. The system is coupled to a small spin reservoir (here 
of size 2). Single spin flips accepted by normal kinetic 
Ising model dynamics are subject to the additional condi- 
tion that the required spin change (±2) can be extracted 
from the reservoir. The dynamics of the system plus 
reservoir are microcanonical in the magnetization, and 
there are none of the questions of energy-transport that 
we have just raised for infinite-range Kawasaki exchange 
p0| . As a result the RG results directly apply and a 
growth-law of L ~ t 1 / 2 is indeed seen ||. 

To check these arguments we simulated critical 
quenches, with (0) = 0, for non-conserved and also 
for globally conserved systems with both Creutz and 
Kawasaki dynamics. For off-critical quenches similar re- 
sults apply. In Fig. [I] we see the expected t 1 / 3 growth 
for T = Kawasaki exchange, and reconfirm i 1 / 2 growth 
for the other cases. To test the rescaling of the kinetic 
prefactor in the Kawasaki algorithm, we show a quench 
to T = 0.9T C and also L/y/p for the quench to T = 0. 
Both of these exhibit the expected i 1 / 2 growth. In Fig. 
||, we see that all of the models have similar spherically- 
averaged scaled correlations. [Anisotropies are too small 
to see for these critical quenches ||.] The correlations for 
the Kawasaki model are slightly different, however this 
seems to be a transient effect. [Strong finite-size effects 
in the T — 0.9T C quench, after the latest time shown, 
limit the duration of the simulation.] 

In conclusion, we see that an infinite-range Kawasaki 
exchange implementation of global conservation laws has 
activated kinetics, with an anomalous L ~ i 1 / 3 growth 
at T = 0. In addition, Kawasaki dynamics leads to long- 
range energy transport. The Creutz algorithm, on the 
other hand, is not activated and has no long-range en- 
ergy transport. It provides a faster and better-controlled 
implementation of global conservation laws, and should 
exhibit the expected i 1 / 2 growth at all temperatures. 
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FIG. 1. Length-scale L(t) for non-conserved (stars), Creutz 
(plus), Kawasaki T — 0.9T C (triangle), and Kawasaki T — 
(circle). The system sizes are 2048 2 , 1024 2 , 512 2 , and 512 2 , 
respectively, with at least 20 samples in each system. We also 
show L/yfp for the T — Kawasaki quench (squares), which 
goes as t 1 ' 2 , as expected. The upper and lower lines show 
and t 1 ' 3 , respectively. 
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FIG. 2. Scaled correlations, divided by the equilibrium 
magnetization squared, from the latest times of simulations 
shown in the previous figure. Solid, dashed, dot-dashed, and 
dotted lines are non-conserved, Creutz T = 0, Kawasaki 
T = 0, and Kawasaki T = 0.9T C , respectively. 
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